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I. INTRODUCTION 


This report presents the results of a ray-trace analysis of several 
glancing-incidence x-ray optical systems. The present study was carried 
out under Contract No. NAS8-32115 by personnel of Montevallo Research 
Associates, Montevallo, Alabama during the period June 7, 1976 to 
A"gust 13, 1976. Participating in the study were Dr. J. William 
Foreman, Jr. and Mr. Joseph M. Cardone. Only Mr. Cardone's time was 
supported by the present contract. Dr. Foreman was supported separately 

uni.',, che NASA/ASEE Summer Faculty Fellowship Program at MSFC. Compu- 

» 

ter time on the Unlvac 1108 at MSFC was also furnished separately to 
Dr. Foreman. 

The object of the present study was threefold. First, following 

up on earlier work*, the vignetting characteristics of the S-056 x-ray 

* % 

telescope were calculated using experimental data to determine mirror 

reflectivities. Second, a small Wolter Type I x-ray telescope intended 

for possible use in the GOES (Geostationary Operational Environmental 

Satellite) program was designed and ray traced. Finally, a ray-trace 

2 

program was developed for a Wolter-Schwarzschild x-ray telescope 
which was designed by members of Dr. A. B. C. Walker's group at 


Stanford University. 


II. V1CNETTINC CHARACTERISTIC; 

S-056 X-RAY TELESCOPE 


The S-056 solar x-ray telescope, a paraboloidwl-li> perboloidal 
Woltcr Type 1 Instrument designed by Mangus and I’t di rv >od i . was 
used as part of the Skylab instrument package to obtain numerous 
photographs of the sun at various soft x-ray wavelengths . In order to 
interpret the resulting photographic data, it is necessary to know how 
the energy throughput of the system varies as a function of off -axis 
viewing angle at various x-ray wavelengths. The reduction it. energy’ 
throughput with off-axis angle (conmonly referred to as vignetting) 
occurs for two distinct reasons: (1) Rays entering the annular aperture 
of the telescope begin to miss either the paraboloidal or the hvper- 
boloidal mirror and are intercepteo by the second stop, which is 
designed to intercept all rays which do not strike both mirrors. The 
second stop is also designed to limit the field of view of t: tilescope 

to approximately * 20 arc-minutes; (2) Owing to variation of the glancing 
angles of incidence at both mirrors v’ith off-axis angle, there is usually 
a small decrease in the overall reflection efficiency of the tc.escope 
with off-axis angle. 

In a previous study , the vignetting characteristics of the S-056 
x-ray telescope were calculated using theoretical mirror r t.iecivity 
data. However, theoretical reflectivity data invariably p. edict higher 
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energy throughput than experimental reflectivity data. It was therefore 
decided to repeat the vignetting calculations using experimental reflec- 
tivity data. Since the S-056 mirrors were made from fused silica, it 
would be desirable to have experimental reflectivity data over a wide 
ran^e of wavelengths for fused silica. However, there is no such data 
of which we are aware. Consequently, it was decided in consultation 

with cognizant NASA-MSFC personnel to use the data for Pyrex glass 

4 

given by Stewardson and Underwood for wavelengths less than 8. 34 A 

and to use the data for F-l glass of Ershov, Brytov, and Lukirskii^ for 

wavelengths of 8. 34 A and larger. 

The reflectivity for a ray which strikes an x-ray telescope mirror 

6 

at a glancing angle 0 is given by the Fresnel equation, 



and 



( 4 ) 



u 


In £q. (4), u c »» the critical angle for total external reflection. 

The critical angle 0 C is given by the relation 

3 C = ( 2 ^) (b) 

where 6= 1 -n and n is the real part of the index of refraction of the 
mirror material. For soft x-rays, n is slightly less than unity for 
all common materials, so that 6 is a small positive quantity (typically 
on the order of 10” 4 ). In Eq. (3), p^ is the linear absorption coefficient 
of the mirror mu.crial and A is the wavelength of the incident x-rays. 

Ershov, Brytov and Lukirskii giv alues of p^ and 6 for F-l 
glass for wavelengths from 8. 34 A to 31. 36 A. These values can be 
substituted directly into Eqs. (1) through (5) to get the reflectivity 71 
for each given wavelength. Stewardson and Underwood, on the other 
hand, give data for Pyrex glass at several wavelengths less than 8. 34 A, 
and they state the experimental values of the parameters 0 c and y rather 
than p^ ~nd 6. By use of Eqs. (3) and (' ), one can easily convert the 
values of 0 c and y into the corresponds g values of p^ and 6: 



z 



where O c in Eq. (7) i a to be expressed in radians. A summary of 
the values of and 6 used at each wavelength is given in Table 1. 

The relative energy contribution from each ray which reaches 
the focal plane of the telescope is R(^p) ’ R(^j_|), where R is the 
Fresnel reflectivity function defined in Eq. (1) and 0 and 0 ^ are the 
glancing angles of incidence at the paraboloidal and hype rboloidal 
mirrors, respectively. The total relative energy reaching the focal 
plane at a given wavelength and at a given off-axis angle is the sum 
of the products R(6p) ' R( °H> for all rays which reach the focal plane. 

The final results of the vignetting computer runs are summarized 
in Table II and plotted in Fig. 1. It will be observed from Fig. 1 that 
vignetting is roughly a linear function of off-axis angle out to approxi- 
mately 20 arc-minutes. Beyond this off-axis angle, the second stop 
begins to come into play, and vignetting becomes more severe with 
increasing off-axis angle. 

At each off-axis angle, 36, 360 rays were entered in the telescope 

aperture. For a variety of reasons, not all of these rays actually 

9 

reached the focal plane. The number of rays actually reaching the 
focal plane at each off-axis angle is summarized in Table III. Since 
each ray entering the aperture is assumed to have an equal relative 
energy weight of unity, the data in Table III can be used to calculate 
the telescope efficiency for any given case. For example, using 



6 


Tables 11 and 111 one finds that the telescope efficiency for an 
off-axis angle of 5. 0 arc-minutes and a wavelength of 31 . 36 A is 
r\ - 17,8 24. 6/ 24, 859 = 0.717 = 71. 7%. All efficiencies estimated in 
this way are probably somewhat too high, since the ray-trace results 
obtained here do not take account of x-ray scattering and absorption 
due to mirror surface roughness, imperfect mirror figuring, possible 
surface contamination, etc. 



7 


Table I. Values of P/ and 6 Used at the Various 
X>Ray Wavelengths as Parameters 
in the Fresnel Reflectivity Equation 







b 


Table II. Result* of S-056 Vignetting 
Computer Runs 


CT-Axis Angle 
(arc-minutes ) 


Wavelength 

(A) 


Relative Energy 
in Spot 


6.16 
6. 62 
6 . 86 

8. 34 

9. 89 
12. 25 
1 3. 34 

14. 56 

15. ‘>7 
17. 59 
19.45 
21. 84 

23. 62 

24. 78 
27.42 
31. 36 


379.0 
1434. 9 
2086. 8 
15535. 1 

17335. 5 

16404. 2 

16784. 6 

16956. 9 
161 15. 5 
16333. 1 
15618. 1 
15405. 5 

15244. 9 

17920. 2 
18142. 1 
19101. 1 


6 . 1 6 
6. 62 
6 . 86 
8. 34 
9.89 
12. 25 
1 3. 34 

14. 56 

15. 97 
17. 59 
19.45 
21. 84 


392.8 
1 332. 5 
2089. 9 

14414. 3 

16147.4 
15293. 5 
15652. 8 
15816. 1 
15033.2 

15237. 6 

14571. 6 

14374. 5 


23. 62 


14223. 5 







Table* 11. (Continued) 


Oft-Axis Angle 

W aveien^th 

Relative Energy 

(arc-minutea) 

(A) 

in Spot 


24. 7H 

l 

16721. 1 

5.0 

27.42 

16929. 1 


31. 36 

17824. 6 


6. 16 

428. 9 

6, 62 

1121.7 

6, 86 

166 3. 3 

8. 34 

l 3263. 0 

9. 89 

15051. 2 

1 2. 25 

14294. 6 

1 3. 34 

14644. 2 

14.56 

14804. 8 

15. 97 

14076. 5 

17. 59 

14272. 1 

19.45 

1 3651. 1 

21. 84 

1 3469. 5 

23. 62 

1 3324. 6 

24. 78 

15671. 2 

27.42 

15869. 2 

31. 36 

1671 1. 1 

6. 16 

383. 5 

6. 62 

867. 1 

6. 86 

1225. 1 

8. 34 

1 1649. 5 

9.89 

13674. 1 

1 2. 25 

13067. 5 

1 3. 34 

13413. 5 

14. 56 

1 3575. 7 

15.97 

12916.4 

1 7. 59 

1 31 03. 7 
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19.45 


1 25 38. 4 







Wavelength 

(A) 


fl'-Axib Angle 
(arc-minutes) 


Relative Energy 
in Snot 






21. 84 

12377. 3 


23 . 62 

12237. 9 

16.0 

24 . 78 

14406. 5 


27.42 

14594.4 


31 . 36 

15373 . 5 

1 

6 . 16 

329.2 


6 . 62 

720 . 3 


6 . 86 

1002. 6 


8 . 34 

10126.9 


9 . 89 

12332. 0 


12 . 25 

11851. 1 


1 3 . 34 

12184. 8 


14 . 56 

12343 . 2 

20 . 0 




15 . 97 

11749. 8 


17 . 59 

11925. 9 


19.45 

11414. 9 


21. 84 

11272.2 


23 . 62 

11140. 8 


24 . 78 

13125 . 0 


27.42 

13300. 2 


31. 36 

14013 . 8 


6 . 16 

218 . 3 


6 . 62 

500.7 


6 . 86 

708 . 3 


8 . 34 

6784.7 

25 . 0 




9.89 

8112 . 0 


12 . 25 

7772 . 6 


13 . 34 

7981 . 6 


14 . 56 

8084. 6 
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FIGURE 1. NORMALIZED RELATIVE ENERGY IN SPOT AS A FUNCTION 
OF WAVELENGTH FOR VARIOUS OFF-AXIS ANGLES. 
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III. DESIGN OF A WOLTER TYPE I X-RAY TELESCOPE 
FOR THE GOES PROGRAM 

It is tentatively planned to place an x-ray imaging system 
aboard a satellite in the GOES program to observe x-ray activity 
on the sun and correlate such activity with global weather patterns 
and communications disturbances. The following data were furnished 
as basic design inputs for the proposed GOES x-ray imaging system: 

(1) The imaging system will be a Wolter Type I 
paraboloidal -hype rboloidal x-ray telescope. 

(2) The maximum glancing angle of incidence on the 
paraboloidal mirror will be 83 arc-minutes (this 
corresponds to a short wavelength cutoff of approxi- 
mately 6. 15 A for a gold-coated mirror). 

(3) The entire telescope assembly must fit in a 
package approximately 30 inches long. 

Using these inputs, we chose the maximum glancing angle of incidence 

on the paraboloid, 0 , to be 83 arc-minutes and the baseline 

r max 

focal length of the system, ^baseline’ to ® inches. We also 

cho«e to make the length of the paraboloidal rr irr jr, Lp, equal to 
2. 0 inches, aiming toward a naraboloidal collecting area of roughly 
4. 0 cm 2 . 

The basic design equations for a Wolter Type I x-ray telescope, 
rewritten slightly from the forme given in Reference 3, are: 


Equation of paraboloidal mirror: 


l 


= f f) 


(B) 


Equation of hypcrboloidal mirror: 


f * * + £ 


(a-O 


- JL = I 


(9) 


V ,t X. 

<k + b » c 


f 

rHiN 

“^&*seuNe ^ 4 ' a MNx^ 

(10) 

f = 

Amiu (***») 

(ID 

* PfA\r4 

= Cf?K\H “ fO/^p 

(12) 

“^Phax 


(13) 

P 

J Pmnx 

= [f (l V* + 

(14) 

c 

- ^&xseut4e / 2. co* (u- 

(15) 


= c |l« (l8 H(ol ) - I'j 

(16) 

v> 

II 

o 

1^ 

(17) 





Collecting area of paraboloid = 


Collecting area of hyperboloid = 





(. 6 ) 

U*) 


* 


Focm. Plans' 


zc 




-c 




/'if* 


The geometrical meanings of the various terms are summarized 

in Fig. 2. Use oi the input parameters 0 max = 83 arc-minutes, 

x, = 25.0 inches, and L_ = 2.0 inches in Eqs. (10) through 

baseline P 

(20) produces the following results: 


p = 0. 0584851950 in. 

a - 12. 5292425975 in. 


b = 0. 8570224335 in. 

c = 12. 5585193203 in. 


p 

pmin 

Z 

pmir. 

Z 

pmax 

P 

pmax 

Z, 

focus 



= 2.4219062894 in. 

- 50. 1 170386406 in. 

= 52. 1 170386406 in. 

= 2.4697309276 in. 

= 25. 1 170386406 in. 

- 0. 735 in 2 = 4. 74 cm 2 
= 1. 99 in 2 = 12. 84 cm 2 . 








la 


Usually the coordinates (p, . , Z, . ) of the rearmost point on 

hmiR nmin 

tlie hyperboloidal mirror in the XZ-planc arc determined by requiring 

that the hyperboloid be just long enough to catch all on-axis rays 

incident on the paraboloidal mirror. If this is the case, then an 

incident ray parallel to the optical axis and striking the paraboloid 

at (p , Z ) will, after reflection from the paraboloid, strike 

pmax pmax r 

the hyperboloid at (p, . , Z . A straightforward application of 

hmin nrrnn 

analytic geometry to this situation gives the result 


where 

c>4 = t — CX “Vavn 2 0 

^ r - Z \>C 


0 s • 


(20 


( 22 ) 

( 23 ) 

( 24 ) 


( 25 ) 


The value of p. can then be determined from the relation 
hmin 




= * 


H(Alr4 


4-aa *2.0 


( 26 ) 
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The values of Z. and p, . found for the GOES telescope from 
hmin , hmin r 

Eqs. (21) and (26) are 

Z L . = 48. 2690250035 in. 

hmin 

p . . = 2. 2873806150 in. 

hmin 

The length of the hyperboloidal mirror, L , is thus 

H 

L t? = Z . - Z, . = 1. 848 in. 

H pmin hmin 

The results of a ray-trace program to compute the rms spot 
radius in the focal plane as a function of off-axis angle are summarized 
in Table IV and Fig. 3. The effective focal length of the telescope 
is 25. 18 inches, and the plate scale in the focal plane is 3. 10 microns/ 
arc-second. 

As a check on our ray-trace results, we have used the empirical 

7 

formula developed by VanSpeybroeck and Chase for the rms spot 
radius in the focal plane as a function of off-axis angle. The 
empirical results are plotted in Fig. 3 for comparison with our exact 
ray-trace results. There is reasonably good agreement between the 
two. 

We have also determined the field curvature at the paraboloid- 

V 

hyperboloid focus, using a ray-tract: program which automatically 
locates the plane of best focus for any given off-axis angle. The 
plane of best focus is defined to be the plane in which the rms spot 





R.M.S. SPOT RADIUS (ARC-SF 



Figure 3. Root-mean-squarc spot radius in the flat focal 
plane of the GOES x-ray telescope as a function 
of off-axis angle. 



'll 


I 

radius is a minimum. Table V and Fig. 4 sumn... -ir.c *he field 

8 

curvature results. The empirical formula for Z, , - Z, , , 

best focus focal plane 

given by VunSpeybroeck and Chase has also been used to calculate 
the field curvature as a check on our ray-trace results. The 
empirical results arc plotted in Fig. 4 and show good agreement with 
the ray-trace results. 

The rms spot radius in the surface of best focus as a function 
of off-axis angle is given in Table VI. Comparison of the results in 
Table VI with those in Table IV shows that the spot size in the surface 
of best focus is considerably smaller than the spot size in the flat 
focal plane at any fixed off-axis angle, as one would expect. The 

i 

two spot sizes are compared graphically in Fig. 5. 
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Table V. Field Curvature Rcauita at the 
Paraboloid-Hyperboloid Focua of the 
GOES Tele a cope 


z 


Olf-Axia Annie 
(arc-minutca) 


8 . 0 
10.0 
12.0 
14 . 0 
16.0 
18. 0 


rz 

Beat Focua Focal Plane 
(millimetcra) 

0. 0000 
0.0038 
0. 0152 
0.0341 
0.0602 
0.C939 
0. 1347 
0. 1810 
0.2346 
0. 2988 


20 . 0 


0. 3662 





Z BEST FOCUS - ^OIAL PLANE (I. ! 1 LIMITERS) 


SOLID CURVE: RAY-TRACE RESULTS 

POINTS S l EMPERICAL RESULTS 

(VANSPEVBROECK A CHASE) 



FIGURE 4. CURVATURE OF FIELD DATA AT THE PARABOLOID-HYPERBOLOID FOCUS. 



ij 


Table VI, Root-Mean-Squarc Spot Radius in the 
Surface of Bc6t Focus of the GOES Telescope 
as a Function of Off-Axis Angle 



20 . 0 


6.15 





OFF-AXIS ANGLE (arc-minuta*) 

FIGURE 5. COMPARISON OF THE SPOT SIZES IN THE FLAT FOCAL PLANE AND IN 
THE SURFACE OF BEST FOCUS FOR THE GOES TELESCOPE. 




27 


IV. IAY-T1ACB ANALYSIS OF A WOLTER -SC.1WARLSCH1LD 
X-RAY TELESCOPE 

The Wo Iter Type I paraboloid -hyperboloid x-rey telescope, which haa 

found auch wldeapread uae In x-mv & • ronony up to the preaent rise, la 

fully corractad for apherlcal aberration but only approximately corrected 

for cooia. It la poaalble to design an aapherlc-aapherlc x-ray teleacope 

which la completely corrected lor both apherlcal aberration and coma. 

An aplanatlc x-ray teleacope of this type la called a Wolter-Rchwariachlld 
2 

x-ray telescope . A nested array of three Woltar-Schwarsachlld (tf-S) 
telescopes has been designed by members of Or. A. B. C. Walker's group at 
Stanford University. Our taak was to develop a ray-trace program for this 
system. 

Physically, a W-S x-ray telescope resembles a standard Wolter Type 1 
teleacope. However, In the W-S teleacope the two mirror surfaces are 
general aapherlca, rather than conic sections aa In the Wolter Type I 
teleacope. The equations for the mirror surfaces In the W-S telescope 
are expressed parametrically in terms of an angle p : 

MIRROR # 1: (Analog of the parab old In the Wolter Type I system) 



where 

s -Van ( p /0 
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nje parsaeUrs Z Q , 

Dt served from Fig* 6 Chat Z la the baa* line focal length of the ijrati 


arid that 


! 


t 


*. +«T‘ 


(*) 


X u Is ths radius at the Intersection of the two mirrors. 

RROR # 2: (Analog of th<- hyperboloid in the Wolter Type 1 system) 


» JL 


CAg 


f ’ 


where 


ia[ i *'* 4- f 4w» 1 (f/»')-A]' + I 

A 4 1 k wk 1 1 J 3 


For the ortermost aet of W-S mirror* in the Stanford telescope, 
Z D * 50.0 :-.ches and X u * 7.0 inches. 


( 30 ) 


(31) 


(32) 


(33) 


(34) 


. (35) 
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The length L 1 of Hlrror # 1 <"•• Fig* 6) was specified by the Stanford 

group to be 10.5 inches, nakinc Zl, nax ■ 60.5 inches. The first Job was to 

determine the corresponding value of ^ max . For this purpose, an iterative 

computer program waa set up to find the value of ^ corresponding to the 

value I, - Z, ■ 60.5 Inches The result was 
1 imax 


P 


ms* 


- 8.3753315183. 


Incidentally, the minimum vl of & is 


Putting 


f. e 


f$* - 7.969610J933 0 . 
max in Eqs . (2/), (33), and (34) gives 


Xi - 7. 35388002 inches 

X 2min - 6.09112390 inches 
Z 2mln - 41.37225999 inches. 


Ihe leng h of Mirror # 2 la thus 

Lo ■ Z_ - Zo . ■ 8.62774001 inches. 

, 2 o *min 

In erder to ray trace > e outermost set of U-S mirrors, we attempted 
to fit Mirrors # l and # 2 with cubic spline functions* 0 in the usual way. 
However, it turns out that Mirror # 1 has a horizontal tangent at the point 
of inter* ection with Mirror # 2, and this causes a singularity in the 
tpjLine fi notion fit at that polr.t. We then had to back up and attempt a 
cpillne fticMon fit which leaves out an infinitesimal region of both 

t j 

mirror? rear their intersection point. Owing to extremely long turn- 
i r >und times on the Univac 1108, we were just able to get a ray-trace 
program t ased on this limited spline function fit completed before our 
time ran put on the present contract. Tin* results of this ray-trace 
prigram «re summarized in Table Vtl, More work is needed to get a better 
spline ft notion fit for the min r < (especially Mirror # l) and to double 
check the results in Table VII. it is also necessary, of course, to 
complete the ray trace for the other two sets of nested mirrors. Unfor- 
tunately, we were not able to complete this work under the present con- 
tract. 









32 


) © 

^rFERENHES 

1. J. W. Foreman, Jr. and J. M. C.udone, Vignetting Character la 1 lea of 
tha S-056 X-Ray Telescope . Final Raport, Contract No. NAS8-31996, 
Montavallo Raaearch Associates, Montevallo, Alabama (April, 1976). 

2. R. C. Chase and L. P. VanSpeybroeck, Wolter-Schwarischild Talaacopaa 
for X-Ray Astronomy . Applied Opt ’’ r a 12.* 1042 (19! f. 

3. J. D. Mangus and J. H. Underwood, Op tical Design of a Glancing* 

Incidence X-Ray Telescope . A t i-d Optica 8, 95 (1969). 

4. E. A. Stewardson and J. H. Underwood, The Reflection of Soft X-Rays 
by Polished Surfaces of Gla and Steel . British Journal of Applied 
Physics 16, 1877 (1965). 

5. 0. A. Ershov, I. A. Brytov, and A. P. Luklrskll, Reflection of X-Rays 
from Certain Substances in the Region from 7 to 44 X . Optics and 
Spectroscopy 22 ., 66 (1967). 

6. The Fresnel equation is writ'-n here in the torra given by R. W. Hendrick, 
Spectral Reflectance of S o! id. for Aluminum K Radiation . Journal of the 
Optical Society of America 47, 165 (1957). Strictly speaking, one 
would expect two Fresnel equations to apply, one for radiation polarized 
perpendicular to the plane of incidence and the other for radiation 
polarized parallel to the plane of incidence. However, for the small 
glancing angles of inciden c used here, the two equations give virtually 
identical results, so that for practical purposes the Incident radiation 
may be considered to be unpolnrlzed, with Kq. (1) giving the reflectivity 
for unpolarized Incident radiation. 

7. L. P. VanSpeybroeck and R. C. Chase, IV.< i ;»n Parameters of Pnraholold- 
llyper hol • ■ c j T« • 1 . ■■■c opes f r X-Ray Astronomy , Applied Optics H^, 440 (1972), 
Eq. (2). The- tir t term in this . qua t ion must be doubled to give the 


OBlGWAtFACTlS 
OP POOR QUALuii 


1 s 


rr.ij sj ot radi in a flat . 
to th< rms sp<" iius in i 
1 . r aent ion . 

6 . Refer* n <• 7 , Eq. (6) . 

S. Ri ftrincc 1, Section II. 

1( . A. i . ler and T. P. V 1, I-'ui k- l i ons : An Alternative Kopro- 

si-nlat to ' ' phcric S ur face ;• 1 i t-vl Optics J_0, 1<>4H (1971). 


' plane. 'lilt! equation as it stands rclrr. 
irved surl.i •’ ol" host locus, as the authors 




